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Exact Treatment of Antenna Current Wave Reflection at the 
End of a Tube-Shaped Cylindrical Antenna 
ERIK HALLEN+ 


INTRODUCTION 


NTENNAS, whose radiating body has a cross 
section which is small in comparison to the 
length of the antenna and to the wave length; 

i.e,, antennas consisting of wires and rods, are to a 
very wide extent treated with the help of a linearized 
integral equation, the invention of the author.!-3 
In this equation the distance between two points on 
the antenna is normally represented by the distance 
between the corresponding points on some central line. 
Only when the distance is small this is not permitted 
and from such regions arises the only term which 
contains the dimension of the cross section, which is a 
parameter mainly consisting of a logarithm. The equa- 
tion therefore has a certain limited degree of accuracy 
which is such that the ratio of the radius of cross section 
to the length of the antenna or to the wavelength is 
neglected compared with unity. The results which can be 
drawn from the linearized integral equation thus also 
should have this limited accuracy which is a normal 
one in electrotechnics in all kinds of devices, where 
wires are involved. Nevertheless much discussion has 
gone on about this accuracy. The only way of finding 
definite numerical answers to this question is to solve 
exactly the antenna integral equations, both the 
linearized one and the exact one, for some special case. 
Nowadays this can be done for a straight cylindrical 
tube-shaped antenna. 


(R. Gans has recently expressed the opinion that 
Hallén’s linearized integral equation should not have 
any exact solution.*> This is a mistake made by Gans 
because he apparently has never seen my original 
papers. What he studies and criticizes is the coarser 
form of the equation given in many American papers 
and books as ‘‘Hallén’s integral equation.’”’ Gans in 
reality criticizes the deviation that is made in those 
papers from my own form, which is not subject to 
any criticism of the kind expressed by Gans.®? 
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It must also be remembered that even if I would 
agree with Gans, to some extent, if he had directed 
his criticism to the proper persons, it is for practical 
purposes not very essential if the coarser form of the 
equation is used in the beginning. Those authors have 
at a later stage in their deductions usually completed 
the linearization so that they in reality use the correct 
form even if it has not been written explicitly. As they 
try to find only an approximate solution it is even for 
this reason unessential if their equation is not quite 
the correct linearized one.) 

Any electric field in free space, periodic with respect 
to time as e, can be expressed with the help of a 
vector potential A as follows: 


E = 5 (grad div A + B?A), (1) 
Jj 
where c is the electromagnetic wave velocity and 
B=, (2) 


c 


a wave constant. On a cylindrical antenna oscillating 
under conditions which involve axial symmetry, the 
current and consequently its vector potential has the 
direction of the cylinder axis. If the ohmic resistance 
is neglected the tangential component of E vanishes 
in the surface of the antenna and from (1) follows the 
nowadays well-known fact!~* that the vector potential, 
and hence the scalar potential, of the antenna field 
is exactly sine-shaped along the surface of a cylindrical 
antenna. This gives the integral equation for the 
antenna current I(x) of a cylindrical antenna, fed by a 
potential jump 2V, in the middle (reference 3, eq. 
35.15): 


=—(- V .e~ %8| =| +A cos fx), (3) 


where 2/ is the length, a the cylinder-radius of the 
antenna, 
r = [(w — 8)? + 4a? sin? $¢]"/?, 

the distance between two surface points, Z, the wave 
resistance of free space (Z, = 377 ohms), and A an 
unknown constant. The minus sign on Y,, which we 
have introduced here, although it was not in the 
previous papers, only indicates that we have turned 
the potential jump so that the outgoing potential wave 
is positive when traveling in the direction of negative x. 
In order that (3) shall be mathematically exact the 
ohmic resistance should be negligible and the cylinder 
should have no end surfaces; i.e., it must be a thin tube. 
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(In practice only the pieces nearest to the cylinder ends 
need to be hollow.) 


The current I(x) can be considered either as a 
standing wave or as a system of traveling waves, one 
outgoing current wave traveling in both directions 
from the feeding point, and a series of current waves 
reflected at the ends of the antenna. The terms of the 
right side of (3) are, besides a constant factor, the 
corresponding terms of the vector potential along the 
antenna. Thus — (Vo/e)e~ #8] >| is the outgoing vector 
potential wave and A/e cos Bx the sum of all the 
reflected vector potential waves. The outgoing scalar 
potential wave is + V.e#* in the direction of negative 
x and — V.e—%#* in the positive direction. The sum of 
all the reflected scalar potential waves is —jA sin Bx. 
Thus the traveling potential waves have constant 
amplitudes (on the surface of the cylindrical antenna), 
whereas the corresponding traveling current waves 
always decrease. 

The outgoing traveling wave is easily found from (3). 
It is the same as the current on an infinite antenna 
without end reflections, hence A = 0 andi] = o. 
Thus the outgoing current wave satisfies the integral 
equation®®, 


[Pea p [ae - 


[For the solution of this equation see (9).] Before we 
proceed we will even write down the linearized equa- 
tions, which correspond to (3) and (4). These are for 
an antenna of finite length? (Eq. 24): 


A(? — x)  1(g)e7 #le-#! — (x) 
a? [a — &| 


= 4a (— V.e~i6l 2! + A cos Bx), 


= $ Voe- ||. (4) 


I(x) log + dé 


(Sa) 
or in a modified form: 


T(x)[log a — 2y — jr — eB 2(B(L — x)) 


—eacenad ay} + f_ MD=IO cates at 


= =(- Ve-i#l2l + A cos Bx), (Sb) 


where the function J is the complex amplitude function 
of the sine and cosine integral as defined (and tabulated) 
earlier by the author.*-!9 From (5b) we immediately 


get the linearized integral equation for the outgoing 
waves, if we put] = © and A = 0.38 
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I(x) [log zm = By = jal 
© (6) — Tle) 
is aa ei 


| x — 

The difference between (6) and (4) is that in (6) 
the radius of the cross section of the antenna is conse- 
quently suppressed except in the first logarithmical 
term. There is the same connection between (5a), (5b) 
and (3). The linearized equations are much easier to 
handle, but that is not the main reason to use them as 
we can now even solve (3) and (6) exactly. The linear- 
izing has been invented in order to make it possible 
to solve antenna problems of much more complicated 
nature. Even if the antenna has a curved instead of 
a straight central line, has variable cross section and a 
cross section which is not circular, integral equations 
can still be set up with only one unknown variable, 
of a type corresponding to (5a). The fact that in the 
general case the vector potential along the antenna is 
not sine-shaped, presents no obstacle. If even the 
ohmic resistance of the antenna is taken into account 
as well as an incoming outer field, this general linearized 
antenna equation can be written :? 

— x2 


=e 
oa ae Pp 
I(x) 


- {| PS - ae 


+f, I(é) af [8 sin B(x — s) 1 — cos (é, s) em tart) 


sit age 8 val 


I(x) igo 


r(&, $) 
+ cos B(x — (+ 3) om as 
- 2° I(s) cos B(x — 9 oe as 


+j . I(s) sin B(x ~ s)zg(s) ds 
An {? , 40 : 
IF E,(s) sin B(x — s) ds+ jz Vosin B| x — hi} 
+ A cos px + B sin px (5c) 
where 2 = 2 log 2//a(x) is a parameter. (5d) 


Here x and é are the lengths of the arc of the central 
line of the antenna reaching from —lI to 1; r(&, x) is 
the straight distance between two points on the central 
line (thus r(x, x) = 0); a(x) is the radius of cross sec- 
tion at the point x, if circular, and the ‘equivalent 
radius” if not; zg is the inner impedance per unit length 
of the antenna wire. E,(x) is the component along the 
antenna central line of an incoming outer electric field, 
if any, and 2V, the driving potential jump [directed as 
in (3)|, if any, applied in the point x = }; ;and A and B 
constants. The integrals with respect to s, which have 
upper limits x, have arbitrary, but constant lower 
limits. All integrands remain finite even when — = 
or § = 


ELECTROMAGNETIC WAVE THEORY SYMPOSIUM. 


Eqs. (4) and (6) for the outgoing current wave both 
can easily be exactly solved. The general solution of 
both includes an undetermined solution to the corre- 
sponding homogeneous equations; i.e., the equations 
with the right side put equal! to zero. In the case of 
the exact equation (4) this eigen-solution (see further 
below) is the wave guide solution for the tube: 


oe (7) 


z. 
I = constant: e*a 


(reference 3, page 410) where the &, are the zeros of 
the Bessel function of order 0. This solution in all 
ordinary antenna cases (A > a), when the wavelength 
is big enough to make 8 < £,,/a, is aperiodic in space: 
frequency below cutoff frequency of the tube. In the 
case of the linearized antenna integral equation (6) 
we have (see further below) a corresponding eigen- 
solution to the homogeneous equation: 


ee fe 
I= constant-et4@V ato, (8) 


where y,; = 1.781072 or log y, = y = 0.577216 is 
Euler’s constant. Eq. (8) represents traveling waves 
in both directions, which in all ordinary antenna cases 
(A > a) would have a very low velocity. These waves 
have no physical significance. 

Thus the general solutions of both (4) and (6) are 
undetermined by terms of the kind indicated by (7) 
and (8). However, if we add the condition that the 
solutions to (4) and (6) should have symmetry with 
respect to the feeding point, so that the two waves, 
going out in both directions, should be equivalent, 
then the solution of (4) as well as of (6) is unique. 
This solution of the exact equation (4) is (see reference 
3, eq. 35.24; reference 8, eq. 8): 


Tee ee oe 
7 2a 
[ eiaz 9 
Ir (a? — B1(aV a — B)K(aVa? — B)°"’ a 
and of the linearized equation (6), 
efae 
4n . j28 ~ da, 
I(x) = ~ 2% yp, fe? 
a , veo J (a? — 6%) lo 08 — 2 _p2) (10) 


where IJ,and K,!! (pp. 78-79) are the Bessel and Hankel 
functions of imaginary argument and I an integration 
path in the complex a-plane which follows the real 
axis from —o to + but avoids a = +8 as if B 
had a small negative imaginary part. If x is changed 
into —x, (9) and (10) remain unaltered. 

Both (9) and (10) are exact solutions respectively 
of (4) and (6), but they still do not represent exactly 
pure antenna waves. Both contain small terms of a 
type similar to (7) and (8) respectively but with the 
difference that x is replaced by | x |; i.e., (9) contains 


11G, N, Watson, ‘A Treatise on the Theory of Bessel Functions,” 
2nd Ed., Cambridge University Press, Cambridge, pp. 1-804; 1944. 
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outgoing waveguide waves (below cutoff frequency) 
and (10) an extra outgoing very slowly traveling wave. 
Those additional parts in (9) and (10) come from the 
poles of the integrand. By changing the path of inte- 
gration we can separate those small terms from the 
rest, which is the true antenna outgoing traveling 
wave, for (9), according to* (Reference 3, eq. 35. 28a) : 


1(e) = — 2 vag (ae)e- 4, (11) 
where 
v(6x) = e8|>| 4 
v 
1 en ialalV/in-w 
‘ J ————— du 
{ ouV1 — wi (abu) Ho) (abu) 
Ll WES av aut, (12) 
wd 1 UVa? — 1H" (agu) Ho!) (apn) ts 


and in the case (10), (reference 3, eq. 35.36): 
e~ i6| 2|‘V 1-u? 


oe du 
az) = e@|-|2.{ [a 
y (Bx) |-| nl ow wt + a log?(4aBy1u)] 


ao ee |[Vea 
aaa 
1UVy? — Afi + — log?(Za8y1u)] 


du (13) 


(I repeat that the minus signs before V, in (3), (4), 
(5), (6), (9), (10), and (11) only indicate that whereas 
the current is counted positive in the direction of 
positive x, the potential jump 2V, is directed so that 
it would drive current in the negative direction; i.e., 
the scalar potential wave is positive for negative x 
and vice versa.) Another way of separating the true 
antenna current from the total current follows from 
(17c) and (17d). In both cases the function (x), 
in spite of the first factor is perfectly aperiodic. 
(This has been diagrammed in reference 3, p. 415). In 
the case of the exact solution (9), and (11, 12) have 
each a distinct physical meaning: (11, 12) is the true 
antenna current on the outside of the antenna tube, 
whereas (9) is the sum of this current and the wave- 
guide current inside the tube. In the linearized case 
there is no such distinction. The new expressions (11, 
12), and (13) still are exact solutions of the integral 
equations (4) and (6) but only on one side of the feeding 
point, either the side of positive x or of negative x; 
for the other side it is necessary to add an eigen-solution 
of the type (7) respective (8) to make the integral 
equation satisfied for all values of x. The expressions 
(9) and (10) are the only symmetrical solutions to 
(4) and (6). We remark that the terms which make 
the difference between (9) and (11, 12) and between 
(10) and (11, 13) contain the smal! factor ap. 

The close connection between (12) and (13) or be- 
tween (9) and (10) is obvious. We have in the de- 
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nominator of the integrands in (12) 


HS) (2)H (2) = Jo?(2) + Y¥.2(2) 


TD = * log? 4 dyig + 22{-—2 4 - (log $y12 
T 


— log? riz)} +o 


and thus we can get (13) from (12) by suppressing the 
powers of a8 but keeping it in the logarithm. This is 
just what we did when we put up the linearized integral 
equation (5a). The numerical agreement between the 
solution (12) of the exact equation (4) and the solution 
(13) of the linearized equation (6) is very high even 
near the feeding point, in all ordinary cases. Only for 
extremely thick antennas is there a difference. It now 
also appears that the consequent complete linearization 
of the integral equation (5a) with suppressing of all 
powers of a (but keeping a in the logarithm), which I 
undertook before!? is the only strict procedure. My 
followers have mostly written a fictitious 


-Ve-p ba 


the distance between a point on the surface and a point 
on the axis, instead of |x — &|, and even used a square 
root expression in the first term of (5a). In doing so 
they have not at all increased the accuracy and in fact 
have made (5a) unsolvable, in a strict mathematical 
sense, because the derivative of the left side then will 
remain continuous at x = 0, but not of the right side. 
Gans, in his above mentioned papers, criticizes this 
inconsistency, which he quite erroneously attributes to 
me. The remedy which he recommends, the changing 
of the right side of (5a) as well, does not seem adequate. 
It is better to return to the original for this integral 
equation method.? (When there is no potential jump, 
Passive antennas, and the antenna is solid, the distance 
r= ~+/(« — £)? + a? may be used, because one can use 
the vanishing electric field along the antenna axis to 
get an integral equation, but then there are always end 
surfaces.) 

From the exact symmetrical expression (9) for the 
outgoing traveling current wave it is very easy to get 
the potentials and field strengths at any point outside 
the antenna and even inside the antenna tube. We only 
have to apply the formula 


ar CO er Birt iat 
x ao f UES dé 
TJ) Oo -0 r 


ae Vo = B*)K o(p a? — B ), p>a@a 


» (14) 
Qeie?T (pV a2 — B)K (a V a? — 6?) B?), p<a, 


where r = V (x — £)? + p? + a® — 2ap cos ¢ is the 
distance between two points, whose cylinder co-ordinates 
are p, 0, x respective a, y, & By putting (9) in the 
integral expression for the vector potential (similar to 
the left side of (3)) we get with the help of (14) the 
vector potential at any point: 


ELECTROMAGNETIC WAVE THEORY SYMPOSIUM 


eee ee KoloVa® = 8) og 
: a KaVat — 8%) 
(15a) 
Biase LMAO), OE Aiea = 8) 
27s oP a = dap <a 
c ws Ta? — Bg? I,(aV of — B?) 
(15b) 
Ave Ve eie, gta: (15c) 


c 


As the scalar potential V = jc/B8 (@A ./dx) we get from 
(15): 


p= Vv: | aca? — Ko pV a? — 8?) 8) 4 a5 
re—B K (av a _ a & 
(16a) 
; jax x V vn? — 2) 
= vif oe 16b 
e ai poa— 8 Tava — 62) Shean panes 
V=- e~iBle], pp =a. (16c) 


x 
"|x| 
As the magnetic field is B, = —0A./dp we also get 


from (15): 
Bf TA NE Daag > 
¢ rVai— 6 Kava — 6) 
(17a) 
Ve # f e@* — Ti(pV a? — 62) 
Bes - TV a? — £? hivaca Oe 
(17b) 
Vo If _ 2? = Ki(aVa eels 
Be) rae al a 
c w/YT Vat — Bp ae 
=a+0 (7c) 
U. Bf a neE=3 - 6), 
B, = Vo? — B 2 
— 8 1,(aVc2 — 8) a 
p=a-—0. (17d) 
The axial electric field component is 
£; = -# (Fhe *4.): 
B\ ax + 8 
E. -%f jus ENO EB 5 ne Ube 
# r K(aV a2 — B*) 
Vo? — 
z= Uf eiaz 2 To(pV at — 8) as BY) 4 da, p<@ (18b) 
Jr I(aVo? — 6) 
E, = 0, p= a. (18c) 
The radial electric field component is Ep = - : 
p 
— => if a Ky (p a — 6?) da, p >a 
a — BK (eV a? — 6) 
(19a) 
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ae 1a 
a aaa if eS I,(aV a — 6?) Soanase 
(19b) 
E=vV ay Babies : Ki(aV a? =O iy, 
: "J TV a? — BP K(aV a — Bp)’ 
p=a+0 = (19c) 
ae ia = Lavo? — 6) = 6) 
aad le a 8° TaVat ~ 8) 
p=a-—Q0O. (19d) 


From (17c) we get the outside current; i.e., the true 
antenna current J, = B, 27a/p.: 


jas ty Bf at OSD, 
6S Oo; 

Za TV a2 — B? K,(aVa? — B) 
(20) 


and from (17d) the inner current (waveguide current) 
I; =—B asa 


At y jab f 

I; Z. uns og) rVai— BT av ab eSBs da. (21) 
The expression (20) is identical with (11, 12) and the 
sum of (20) and (21) is identical with (9). In (17d), 
(19d), and (21d) the integral can easily be carried 
out and gives a series of wave guide waves below cutoff 
frequency (see reference 3, eq. 35.25). The character 
of the potential jump 2 V, in the feeding point and the 
direction we have given it in this investigation is clearly 
seen in (16c). In the following we will also give some 
other forms and asymptotic values for some of the 
quantities (15)—(21). 


giz 


I(aV a? — 6?) 


INTEGRAL EQUATION FOR REFLECTED ANTENNA 
CURRENT WAVES 


After this survey of the outgoing current wave we 
turn to the main object, the reflected waves. The first 
investigation 1 made on them! was based on the prin- 
ciple that a standing wave expression was first obtained 
and then this standing current wave was dissolved into 
a system of traveling waves. Thus the ‘expression 
(reference 9, eq. 55) for the first, second and following 
traveling waves was obtained. However, there have as 
yet only existed approximate expressions for the stand- 
ing waves, consisting of a series of a limited number 
of terms, obtained by iteration from the integral equa- 
tion (5a). Thus the traveling waves will also be known 
through series with a limited number of terms. This 
method is laborious, as already the third term in the 
standing current wave expression (reference 9, eq. 30) 
is extremely complicated. A new direct method of 
finding the traveling waves has been developed later.’ 
It consists of splitting up (5a) into a series of integral 
equations, one for each traveling wave. This can be 
done by studying the periodicity of both sides of (5a) 
with respect to x as well as with respect to /. Thus the 
series of integral equations (reference 3, eqs. 35.61, 
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35.62a, b, c) for the outgoing wave, for the first, second, 
etc., reflected waves, is obtained. From these equations 
it is very easy to find the result (55) in reference 9, 
and even increase the number of terms; i.e., the ac- 
curacy, in the series solutions (reference 3, eqs. 35.65, 
35.66). 


These new integral equations also invite us to try to 
solve the antenna problem exactly in the form of in- 
tegrals without any abbreviated series. In fact, since 
Levine and Schwinger?” treated sound waves in a limited 
unflanged pipe with the help of Hopf’s and Wiener’s 
solution to a certain integral equation, the integral ex- 
tending from 0 to ©, it has been evident that the 
problem of the cylindrical antenna (with tube-shaped 
ends) can also be solved mathematically exactly. In 
my book? I have also given the exact solution of the 
problem of the first reflected wave, as it is derived from 
the linearized antenna integral equation (5a), when the 
distance / from the end to the feeding point tends to 
infinity. [t is also quite certain that a finite / makes no 
real extra difficulties and that the following waves can 
be exactly determined as well, although formulas be- 
come more involved in these cases. It is our object in 
this paper to make a corresponding investigation based 
on the exact integral equation (3) for a tube-shaped 
cylindrical antenna; i.e., to obtain the exact expression 
for the first reflected wave, and to compare the result 
with the corresponding result from the linearized equa- 
tion. In doing so we will get an exact knowledge of the 
particular role played by the antenna ends. Of course 
the tube-shaped ends are only one of several possibilities, 
but for the role of the ends the tube end can serve asa 
type. As I expected—I do not answer for others—the 
difference between the results from the exact equation 
(3) and the linearized equation (5a) consists in terms 
which have the small factor a8. For thin antennas in 
fact the linearized equation is the most general, it is 
independent of the particular form of the ends, and what 
is common to all antennas with differently shaped end 
surfaces is just what the linearized equation gives. 
But it is very important to know, numerically, for a 
special case how far the agreement between the exact 
solution and the solution of the linearized equation 
extends. 


Two features can be expected to be of special interest. 
When the reflected traveling current wave has traveled 
a considerable distance from the end it is obviously of 
very little concern whether we have used the exact or the 
linearized equation. But in the beginning and especially 
at the very end-point of the antenna the difference must 
be expected to be the greatest. Therefore we will deter- 
mine the exact end admittance of the reflected wave as 
it is derived from (3) and compare it with the end 
admittance derived from (5) which has already been 
determined.’ Further, the open tube-end will carry a 


12H, Levine and J. Schwinger, ' ‘On the radiation of sound from 
an unflanged circular pipe,” Phys. Rev., vol. 73, pp. 383-406; 1948. 
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small internal current, which analog to the inner 
current of the outgoing wave can be expected to be 
aperiodical in space and only represent an end-charging 
current. This mainly capacitive current branches off 
from the reflected wave current at the end of the antenna. 
We will determine it numerically at the end-point. This 
current has no correspondence in the linearized solution, 
where there is, however, just as was the case with the 
outgoing wave, a fictitious small, slow, extra wave current 
without physical significance. Besides this, the exact 
shape of the field round the antenna end will be of 
interest. 

For the deduction from (5a) of the complete system 
of integral equations, one for each traveling current 
wave, see reference 3, pages 421-425, 433-434. Although 
the main feature of the reasoning is based only on the 
true antenna waves, it can be repeated in exactly the 
same way even if we start with the exact standing wave 
equation (3), in spite of the fact that its general solution 
contains even the aperiodic tube waves. In this paper 
we will, however, limit the investigation to the first 
reflected wave, and, for simplicity, assume the antenna 
length 2/ to be very great. The influence of J on the field 
and current conditions around the end apparently is 
very small anyway. Then we do not need the whole 
system of integral equations and can obtain the needed 
one in a simple way. 

We introduce in our equations the distance z along 
the antenna from the end and thus have x = z —1, 
where z = 0, when « = —/. We usually need not extend 
2 to greater values than z < / and (for this part of the 
antenna) have | x | = 1 — g. The incoming current wave 
according to (9) is in our new co-ordinate: 


- iyi 


Zo 2a 


L(g) = 
ei) 
/ = da. (22) 
r (a? — 6*)I.(aV a? — B)K.(aV a? — 6?) 


If 7? increases infinitely there will be no reflected wave 
coming from the farther end z = 21, but only, besides 
the incoming current I(z) and corresponding potential 
wave V,e~#8!'-4l, one reflected current wave iz(z) and a 
corresponding potential wave Vyie—%8*, Thus we get 
from (3), when/— o theintegral equation 


[ ne) +agnlar Af” 


= > (—Ve7i8l 21 4. Vye~ #82), (23) 
E+ landr = V(z — £)? + 4a? sin? do. 


Now, if we fictively define Z(¢) as before, even for 
negative ¢, we have, according to (4), 
= Ve~ 781-41, 


eo er 4 
do yard fo ——dg=-F 


Subtracting this equation from (23) we get the integral 
equation for the first reflected current wave 


Ipr 
<— de 


where ¢ = 
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2s om iBr 


Os, 1 4a : 
i a(t) df — f de = \— Vye~ #8, g > 0 (24) 
2 2r Jo Y Zo 
g(z), 2 <0, 


where 7(z) is, for z > 0, the unknown reflected current 
wave, and, by arbitrary definition, for z < 0: 


uz) = —I(z). (25) 


g(z) is an unknown function which we need not deter- 
mine. V: is a potential amplitude, which has to be deter- 
mined. Eq. (24) is valid whether the antenna length is 
infinite or finite. Now we will use the simplification, 
which comes from making / = «. The only consequence 
of this is that the incoming current (22) will be simplified. 
In fact, the expression (22) for the incoming current 
wave approaches, for big J, 


I(2) = — z Vie Ate-OQ,_ 1 (26) 
where 
2(L — it 
G4 = og UO GE a teat = aly 1093 
a’B 2 


In (26) terms of the order of magnitude 2;-% have 
been neglected. Neglecting even 2 beside / we get a 
constant amplitude of the incoming current wave by 
substituting 2; for Qi. , where 2, = log 21/a*B — y — 
ju/2. The corresponding asymptotic expressions for 
the incoming wave of the vector potential and the 
scalar potential as derived from (15a), (16a) are (with x 
replaced by z — 1): 


A; = — Vo eff) — Vo Boe 2 city Pm p >a 
c 
(28) 
V = V,e%8(-9 4. V,2 log © ef8(-99, 1, p > a, (29) 
re) 


where / is very big in comparison to 2, p, a. For p < a 
we simply get from (15b) and (16b) the asymptotic 
values: 


— Ve piste 


c 
V = Verh, (31) 


and in this case the approach is closer; the neglected 
terms are small essentially as e~”/*. 

The first terms in (28), (29) together with (30), (31) 
form a potential ‘‘field” of constant amplitude and 
constant direction. It has no physical content because 
the field strengths B and £ as derived from it are zero 
everywhere. The electromagnetic potentials of a radia- 
tion field are more purely mathematical and less physical 
than in quasi-stationary cases. It is therefore more 
adequate to use the incoming current wave amplitude 
as a measure of the incoming wave than the potential 
[which we keep only in the unimportant first terms of 
(28)-(31)]. We put the constant 


T= 7 Ve i891 
Zz 


o 


A, = (30) 


(32) 


ELECTROMAGNETIC WAVE THEORY SYMPOSIUM. 


and have for the incoming wave the asymptotic for- 
mulas 


I(z) = —Ieib= (33) 
A,= _ Vo e7Blz-F) Hel toe? este p>a (34a) 
Cc 27 p 
_ Vo - 
A,= ef8lz-), pn Ca (34b) 
c 
V = V,e78(2-) + <= og © 7 p>a (35a) 
V= a es p<a (35b) 
B, = bol eiB: p>a (36a) 
2p 
B,=0, p<a. (36b) 


The expressions (34)-(36) have an independent physical 
sense only for z > 0. Forz < 0;i.e., outside the antenna 
end, one cannot reasonably make a difference between 
incoming and reflected wave (or this difference would be 
artificial). In that case (34)—(36) represent a first part 
of a total wave (see further below). 


SOLUTION OF THE INTEGRAL EQUATION FOR THE 
REFLECTED CURRENT WAVE 


We now introduce (33) in (25) and (24) and turn to 
the exact solution of the integral equation (24). The 
procedure for solving (24) is exactly the same as that 
which I have carried out before for the corresponding 
linearized equation.® The only difference is that we have 
a different kernel, which now is the exact one for a tube- 
shaped antenna. The mathematical method is in its 
main feature that of Hopf and Wiener!’ and it was first 
used on a physical tube problem (an accoustical one) 
by Levine and Schwinger.!* Our procedure is, however, 
entirely different from that of the latter. One of the 
advantages we get is that we never need a study of the 
field to get boundary conditions, but that these are 
automatically filled and all constants known. [See the 
very simple deduction of the formulas (48)-(50) below]. 
In fact, we do not need the formulas for the potentials 
and fields (15)-(19), (34)-(36), etc. for solving the 
problem; however, we get them, as we have seen, very 
easily and they are of interest in themselves. 

In order to facilitate our computation we temporarily 
assume that 8 has a small negative imaginary term 
8 = 8: — je so that 8 is situated in the lower complex 
semi-plane. We dissolve the kernel of the integral equa- 
tion (24) into a Fourier integral (reference 3, eq. 35.22): 


1 {?" e— 78r 
fle : 
1 7.2) 
= on | ue CMa a? = B)K (a a? — 6") da. 
Even i(¢) a expressed as a Fourier integral 
: cr Poca 
uy) = aa et@$[Fi(a) + F_(a)] da, (37) 
Qe J —o 


E.C, Titchmarsh, ‘Fourier Integrals,’’ Oxford University 
Press, Oxford, England, pp. 1-390, especially p. 339; 1937. 
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where 

co. 
F,(@) -/ t(uje—@* du, 

oO 

0 
F_(a) -/ t(u)e— 4" du- F,(a) 

—co 


is unknown but regular in the whole /ower complex 
a-semi-plane. When a approaches —j© we get the 
limit value of F4(a): 


4(0) 


F(a) — ico) e-ie4e dy = os (38) 


F_(a) is known because for negative u, according to 
(25), (33), we have z(u) = Fe##“, and thus 


® I 
= rf ei(B—a) u du =———-» 
-% j(B — a) 


which is a regular function of a in the whole upper com- 
plex a-semi-plane. We finally dissolve even the right 
side of (24) and find it to be 


x ia e#[G.(a) + G_(a)] da, 


where G,(a@) is known: 

An An 
G = 
a= ft Ze Ge 


eatin Vie- iB tg— jau dy _ 
which is a regular function of a in the whole lower 
complex a-semi-plane. G_(a) is unknown: 


0 
Ga) = I. g(uje—te" du 


but regular in the upper complex a-semi-plane. It 
disappears at infinity essentially as 1/—ja in this semi- 
plane. 

Introducing the Fourier integrals in (24) we get 


i Wak fo entry + Ral di 


wt 
2a 


F_(a) (39) 


, (40) 


tone ae (av a? — B? — B)K, (aV a? — 8) B?) da 


a jae a a)| da 
is [Gi(a) + G_(a)] da. 


Using again Fourier's integral theorem in the op- 
posite direction the integrations in ¢ and in y give us 
back [F,(a) + F_(a)] and our integral equation attains 
the final form 


it a eta} 27 (aV we B2)K (av a? — f?) 
rj] —-2 


‘[Fy(a) + F(a)] — Gy(a) — G(a)} da = 0. (41) 
Following the Hopf-Wiener method, we write 
21,(aV a? — B)K(aVar — p) = 24 (42) 
x(a) 
where 
oe Lf. ae 
gi(e) = exp ig Ree er Saar 


- log [22.(a/ 5? — B)K(aV/e? — 6] (43) 
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(a) = exp : 
e Qaj J ~w4tjef — a 


log [22 (at? — B)K (a/e? — B)]. (44) 
Here 5 is a positive number so small that 6 < «¢, ie., 
both the integration paths in (43), (44) go between 
B = 8, — jeand —8 = —£, + je. Then g(a) is regular 
and different from zero in the whole upper complex 
e-semi-plane, g2(a) in the lower semi-plane, and both 
these regions can be extended somewhat to the other 
side of the real axis provided that the border line remains 
within the strip +6. We even find 


e1(— a) g2(a) = 1, (45) 
Eq. (41) now takes the form 
1 {* —., ela) fo? — B? 
aes jaz FAT, Fila F(a 
BT om PEF re) + FC 
i IT tay a: G-(a)) i=. (46) 
pila) 


Here gi(e)/(a? — B?), (a? — B)/go(a), (a? — B?)/ei(a), 
F_(a) and Gi(q@) all remain analytical within the strip 
+6. If we anticipate, as is natural, that when 6 is 
slightly complex (8 = 81 —~ je), i(z) decreases at least 
as e-“*(z > 0), and g(z) at least as e“(z < 0), then 
even F,(a) and G_(a) will remain analytic within the 
strip. The integration path in (46) thus can be arbitrary 
if it remains within the strip. The necessary condition 
for this is that the bracket expression vanishes and 
therefore with respect to (39), (40) 


POs | rag ey 
goa) g2(a) 
.a—B4rVi , & — # 
a dass G_(a@). 47 
oO me ae. ee SEP 


The left side is analytic in the lower semi-plane, the 
right side in the upper semi-plane; within the strip 
both sides are equal. Consequently they represent the 
same analytic function. As F, and G_ are bounded, 
P(a) must be a polynomial of no higher than the second 
degree, and as Fy and G_ vanish at infinity it cannot 
even be of more than the first degree. We may thus put 


P(a) = Ci + Cole — 8), 
where C, and C2 are two constants. 
Putting a = 8 we get, from the left side of (47), 
. 28 
Ci =f i. 
err) 


If we divide (47) by (@ — 8) and let @ go towards 
—je, we get with respect to (38), 


ats. : a+, IF — i(0) 
= SFP jg) 43-2 Ppa jp 
joape(a) (a — 8) ¢2(a) g2(a) 

and thus according to (25), (33): 
C, — 0. 
Hence P is a constant: 
. 28 
P= I. (48) 
2 ¢2(B) 
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We determine even V; by using the right side in (47) 

where we put a = —8, which gives 
= j 28 4rV, 

¢1(~—8B) Zo 
or with the help of (45) and (48), 
= So 
Amgs"(B) 
We can even introduce the end admittance of the reflected 
wave (when the incoming wave comes from infinity): 


(49) 


t 


judo Ames 
Yo Vi Z. 2 (8). 
(See Fig. 1 for numerical values.) 
The Fourier transform of the unknown function 7(z) 
is now known, for from (47), (48) we get 


. 26 g(a) 
Fila F(a) =j—TI A 
+ (a) + (a) J on(B) a? — B 


From (37) we now find the solution of the problem. 
The reflected current wave is 


(50) 


($1) 


i(2) = 15 - ae ies (52a) 
or, with respect to (42), 
jose 18 

2792(8) 


(i a e187 1(a) 

: —— —=———- da 
-2 (a? — 8%) ,(av/a? — B’)Ko(av/ a? — B”) 
The factor before the integral in (52a) can also be 
written V; j48y2(8)/Z. . The properties of ge2(a), z(z) 
and g(z) will be studied in the following. 

We state that all the equations (37)—(40), (45)-(52a) 
are exactly the same as when the equation was linearized. 
Reference 3, eqs. (35.52)-(35.57). Thus (52a) and (52b) 
even represent the reflected current wave, as derived 
from the linearized equation, but with gi and g: defined 
as in (43), (44) with the difference that the bracket 
[27.K | is replaced by a logarithmic expression, which is 
the limit value of the Bessel function expression if 
powers of a@ are suppressed. Thus there is exactly 
the same relation between our solution to the reflected 
wave integral equation (24) and the corresponding one 
of the linearized equation as there is between the two 
exact solutions (9) and (10) of the exact integral equa- 
tion (4) and the linearized one (6) for the outgoing 
wave. It is however our task to carry out the comparison 
completely numerically and find out how thick the 
antenna may be if the linearization shall still be per- 
mitted. 

It is easy to check the result (52a). For  < 0 we can 
add an infinitely big circle in the lower complex e-semi- 
plane to the integration path. As ge(a@) is regular in the 
lower half-plane, the integrand has no other singularity 
within the contour than the pole a = B = 61 ~ je 
and we get for negative z the value 


ie) = 128 (api 92) gies = Joins 
Magy ag 


(52b) 


z2< 0, 
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9.010; 
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0,008) 


0.006 


0.004) 


0,002 


40 


Fig. 1—Real and imaginary part of end admittance of tube- 
shaped antenna. Solid black line indicates exact value: (67) and 


in agreement with (25), (33). This is the fictive current 
which we have introduced to compensate the fictive 
continuation of the incoming current wave: 


t(2) + [(z) = 0, ¢ <0. 


For z > 0 we introduce (52b) in the left side of (24) 
and find this side of oe to be, with respect to (14), 


oe da 
nes a =f 
and (for positive 2) we close the integration path in the 


upper semi-plane, where ¢:(a@) is regular. The only 
singularity within the contour now is the pole a = 


(53) 


—B8 = —B, + je and (53), on account of (45), (49), 

becomes 

p28) e~iBe = a e-iBz = 4 Vie 787, ¢ > 0, 
¢2(B) g2°(B) Zo 


in accordance with the right side of (24), which equation 
thus is fulfilled. 

For negative values of zg the expression (53) represents 
the hitherto unknown function g(z) of (24). 

The reflected current wave (52), with z > 0, consists 
[as was the case with the outgoing current (9)| of two 
parts: an outer antenna wave 7, and an inner wave 
guide wave 1; , below cutoff frequency, which penetrates 
into the tube from the end but which decreases ex- 
ponentially. For thick antennas the latter is not quite 
unessential as it represents a current which branches 


(72). Lighter line indicates linearized value: (68) and (69). Broken 
line indicates approximate series (71). 


off from’ the main antenna current wave. It is mainly 
a charge current, and if the antenna has, for instance, 
flat ends instead of tube-shaped ones, it corresponds 
to the small current which charges and discharges the 
end surface. We can easily split (52b) into these two 
currents. We observe that 


1 
I.(aV a — B)K(aV a2 — 6?) 


= 2 2 qi x] 

° i E Tk, aes (54) 

and get 

(2) = a2) + H(2), 

where 

Fee pee etp(a) | Ki(aVe? — 6), 

€ 2792 () -o Va? — Bp K ves 

z2>0 (55) 

$e Sr = jag e720, (a) we (avo? — B?) 2 _ B2) 

t 2792(B) o Woe — Bg? c= Be T (av a 2 “— B) , 
z2>0. (56a) 


With residue calculus we find the latter expression to be 
the following series: 


Se 


i(z) = JI fo8 


VE 9.20262 m? 028 22lay, (j/a V — a8?) 
¢2(8) “T 


Vint — a 


? 


(56b) 
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where the é,, are the zeros of the Bessel function J,(£). 
(56b) is a series of waveguide waves (below cutoff 
frequency) to which we will return later. For a corre- 
sponding series expression for the wave guide waves 
of the outgoing antenna see (reference 3, eq. 35.25). 
The linearized equation does not give these waves (56b). 


THE ELECTROMAGNETIC FIELD AROUND THE 
ANTENNA END 


We are now in the position that we can give an exact 
description in every detail of the field around the 
antenna end (if this is tube-shaped). The way which is 
physically most straight forward would be to use the 
true currents, the incoming current (33) and the reflected 
current (52), for s > 0 and express the vector potential 
from them in the form of an integral extending from 

= 0toz = & [cf. (23)]. But it is easier, and equally 
correct, to extend both currents to z = — ©, with the 
provision (25) that ¢(z) + I(z) = 0 for z < 0. The 
field from the (extended) incoming current J(z) we 
have already in (34a), (34b), (35a), (35b), (36a), (36b). 
The vector potential from the (extended) reflected 
current ae is at an arbitrary point p, 0, z: 


—Hoy afi of” e738 iB Par [* Sah 
Ar are -0 a? — p? 


where 7 = [(2 — 6)? + p? + ee — 2ap cos ¢]!/2. Using 
(14), (42) this gives 


_ 7 edB cpa) KalpV a? — 6%) 4 
oT eee) | 0 ak — 6 KlaVar = Bi) 
p >a _ (57a) 
oe ee e*9,(a) Io(pV a? — #) 
4m’ p2(B) J -20 a? — BOT (aV a? — i) ° 


p <a. (57b) 
For z > 0 this represents the reflected vector potential 
wave, whereas the incoming vector potential wave is 
given by (34a), (34b). For z < 0, i.e., outside the antenna 
end, we can no longer distinguish between incoming and 
reflected waves but we get the fotal vector potential 
wave by adding (34a), (34b). In this region (g < 0) 
it is better to go back to the function ge instead of ¢1. 
We thus get the total vector potential outside the antenna 
end: 


A, = — Ve eit — I? log © eit: 
c 2a p 
HofB_ e®*ea(a) 
277 02(8) J -2 a? — f? 
- To(aV a — B)K.(pV a? — B*)da, p>a (58a) 
, ee Voy iete-n +1 HoJB © e782 p9(c) 
¢ 2r*ys(8) J —» a? — B? 
-K,(aV a? — B)I.(pVa? — B*)da, p <a. (58b) 


The expressions (58a) and (58b) are in themselves valid 
even for z > 0, although we prefer to use yg; instead of 
v2 when zg > 0. The reason why we always make this 
distinction, using ¢: when z is positive but ge when it 
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is negative, is that we can deform the path of integration 
in the upper half-plane in the first case, and in the lower 
one in the second. It can be shown [cf. (64), (65a)] from 
the integrals in (55) and (57a) that they represent, for 
z > 0, waves traveling in the direction of positive z 
(i.e., have a periodic factor e— 48), whereas (52a) for 
z < 0 represents waves traveling in the negative direc- 
tion (periodic factor e787). 

From the vector potential we get the magnetic field 


B, = —0A./0p and thus, using (58a,) (58b) the total 
magnetic field outside the antenna end (zg < 0) becomes 
B= —1 #2 ei: 4 7H [ eitton(a) 

e 2m p 2m?yo(B) J -m Vg? — Bt 

- I(aV a? — B)Ki(pV 0? — B?)da, p>a (59a) 
paw yp _hed® [? eBtexle) 

° In2ga(B) J -~ Va? — B 

K(aV a? — 6?) — 6) Iy(pV a? — 6?) — B)da, p<a. (59b) 


These expressions for the total magnetic field are still 
valid even for zg > 0, but in this case it is better to 
write them [using (42)] as: 


B, Sa I Bo ei8z + J Hoff 
27p 42°”(B) 

% eltpi(a) Ki(pVa? — 8%) 

-2 Vo? — 6 K,(aV a — 6?) =a p>a (59) 
ae of 2 ety (a) LV = 8), 
¢ 4n°02(8) J -0 Vo? — B? Io(aVo2 — B) 

p <a. (59d) 


In the term outside the integral we recognize the field 
of the incoming wave (36). If we put p = a, we find 
from (59d) the limit value of —27a/u, B,(a) to be 
equal to the inner current (56a), and from (59c) the 
outer limit value of + 27a/y. B,(a) to be equal to 
the sum of the incoming current (33) and the reflected 
outer current (55). Our interpretation of the physical 
significance of the two terms (55) and (56a) of the 
reflected current is thus justified. 

It is easy to verify that outside the antenna end the 
magnetic field is continuous at p = a. If we subtract 
(59b) from (59a) and put p = a we find from (reference 
11, p. 80, eq. 20) the result zero. However, it is a lack of 
elegance not to have a common formula for p < a as 
well as for p > a outside the antenna end, because the 
limit p = a@ between the two regions has no physical 
significance. We find a common formula by deforming 
the integration path in (59a) and (59b) into a double 
line in the lower complex a-semi-plane from a = 8B 
to a = B — jx. If we further substitute Va? — g? = 
—jBu we get, for z < 0, the common formula: 


ae HojB i eb u-19,(— 56 Vu — 1) 1) 
. 2rex(8)Jo  Vyr—1 


-J (apu)Ji(pBu) du, (59e) 
valid for all p but only for negative z. The integration 


B 
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path should avoid the singular point « = 1 in the upper 
complex semi-plane. 

We finally give even the components of the electric 
field which are equally easy to deduce. The fofal electric 


field is 


Be fe K eV a? — 6) 
B= — 12 I etare(a) RB 
47? p2(B) — K,(aV a? — B?) a 
p >a _ (60a) 
mn ; I.(pV¥ a — ?) 
E, = — 12 f Jaz ee ee 
aise "LV aoe” 
p <a (60b) 
Boe Zor a Z oj °  aei*7p;(@) 
s 2rp 4? go(B) J -0 Va? — B? 
Ki(pV o — 
Ki(pV a? — 6) ee MIS) 


Kav o? — B?) cs 
BE, = y—Z4_ is aelroi(a) Ti(pV a? — 6) 4 
* Ant po(8) J -0 Vat — BT (aV a2 a 


p<a,. 


(61b) 


The first term of (61a) is the incoming wave (which 
has no g-component). The equations (60)—(61) are 
valid both for positive and negative values of 2, but for 
negative z it is better to write the common formulas 
for the total field outside the antenna end: 


Z 8 © yebe Ve 1 
2rex(6) Jo Vy2— 1 
-g2(—jBV/u? = 1)Fo(aBu)J o( pBu) du 
= 2b . ebevie—1 
27¢2(8) 
‘g2(—jBV u? — 1)J.(aBu)Ji(pBu) du, (61c) 


E,= —I 


(60c) 


where the integration path avoids the point u = 1 in 
the upper complex w-semi-plane. The formulas (60c), 
(61c) are valid for all p but only for negative z. 

Eqs. (52)-(61) give a complete mathematical solu- 
tion to the problem under investigation. 


INVESTIGATION OF THE ANTENNA CURRENT 


The inner current (56b) is aperiodic and decreases 
exponentially. Its initial value at the antenna end has 
most interest 


0 = 7 108 or C/V Em - av B*) 62 
oe ae Vint — op So 


to which we will return later. It constitutes a branch-off 
mainly capacitive current at the end. With the help 
of (49) we may define an inner end admittance 


y= OO = Fi gil?/oW En? — 078? — 2°?) 
- jaBer 0) > Er = ae 


The outer reflected current (55) is periodic and repre- 
sents a wave: By changing the integration path through 


(63) 
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the upper complex a-semi-plane we can write it in the 
form 


tele) = _ ——~ Wn (Ba)e~ 18° = z Vihoo (B2)e#0*, (64) 
where 
Won (82) = — 28 e2(8) 
—0 eile+B)2y,(a) da 
i — — (65a) 
~s (a? — B2)K,(aV a?— B)K.(aV a? — B? Liv) 


(here \/a? — 8? is positive along the integration path). 
By a suitable change of variable this also becomes 


Vo (Be) = n(6) 2 
if eb i-V 1D (84/u? — 1) 
0 


uv/w — 1H (apu) H,') (aBu) 
where z shall avoid the point u = 1 in the upper com- 
plex u-semi-plane. According to (65a) the amplitude 
function W.(6z) is an aperiodic function of Bz. The 
index c is used in order to remind us of the fact that our 
reflected wave arises from an incoming current wave 
which has come from infinity. If it comes from a finite 
distance / the problem can also be solved but gives a 
more complicated formula. A finite 7, however, influences 
the reflected wave rather slightly. It should be noted 
that (65a), (65b) differ only by the two factors g2(f) 
and g, from the corresponding formulas for the outgoing 
wave (reference 3, pages 411-412). 

The exact numerical evaluation of (82) from (65a) 
or (65b) for all values of z is unhappily very complicated. 
One can fall back on series expansions (see further below) 
which are especially good when z is not too close to 0. 
For z = 0 we have however the value of the reflected 
current according to (49): 


i(0) = 4,0) + 4,(0) = I= > (BV, 


du, (65b) 


(66) 


and according to (50) the total end admittance for the 
reflected wave 
_ 10) — 4r 
Vi Zs 
where 2(8) is simpler to determine numerically. That 
has been done in this investigation for a series of values 
of a8, which give the exact values of the end admittance 
of the straight cylindrical tube-shaped antenna, shown 
as a diagram in Fig. 1 and Fig. 2. Instead of a8 we have 
plotted the ratio \/a between wavelength in free space 
and radius of cross section. We have the simple relation 
A/a = 2r/aB. Even Y; has been numerically computed 
for a few values and as soon as sufficient computing aid 
is available to me I will publish a complete diagrati. 
Thus, Yo , Y, and Y; are all known. For details of thé 
numerical determination of ge see below. 

Before we proceed further we will even mention the 
corresponding result as derived from the linearized 
integral equation. In that case there is no inner current 
4;. In fact, all the equations (56b), (62), (63) above 


2'(B), (67) 
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Fig. 2—Rea] and imaginary part of end admittance of tube- 
ON antenna, exact value. Dots indicate inner end admittance 
63). 


contain the small factor a8, which is supposedly negli- 
gible at the linearization. Thus with linearization we 
have (reference 3, eq. 35.57): 


Yoo (0) = ¢2°(8) 


and 
nn -F g2%(8) (linearized equation), (68) 
where in this case 
rte oi 4 
= log log ——__--— 
p2(@) a ae oe og Be Bye 
(linearized equation). (69) 


The function (68) has formerly been computed (reference 
3, page 432, Fig. 412) and it is introduced in our Fig. 1 
here also for comparison with the exact value for the 
tube-shaped antenna. 

For arbitrary values of 2 the amplitude function 
V oo (82), as it follows from the linearized formula corre- 
sponding to (65b) or (65a) [or (64) and (55)], can be 
expanded into a rapidly decreasing series: 


oo (83) ‘mane 0,71 + aa[ + jo. J? — =] 
6 282 


+ Q-F—3p ~2po1 — Jou 4 731 4 31% — fe) 
+ 1,(/! — 1) ~ 2S3]o8. + +++ (linearized equation), 


(70) 


where 2, = —2 log a6 ~ 2y — ja + 1,(262) + 1(2Bz) 
and the different functions / are the amplitude functions 
of the iterated sine and cosine integrals of different 
orders as defined in (Hallen,'* page 4; reference 9, 
pages 4-5; reference 3, page 435; reference 10, pages 
4-6). They are tabulated in (reference 10, Tables IV-VJ). 
The parameter Q, increases with increasing z and the 
expressions within the brackets in (70) decrease, with the 
exception of the constant, essentially as 1/8z. Thus the 
accuracy of (70) increases with g, and it is a very good 
expression. An exception is the very neighborhood of 
the starting point (antenna end) z = 0, where the series 
gives 


2 
vo (0) = 2-1 + 7 Q5* + 450-4 + +: (71) 
ao 
where S; = > 1/n® = 1.2020567 andQ, = — 2 log ag — 
1 
2y — ja. Even this approximate value (the end ad- 


mittance derived from it) is entered in Fig. 1 for com- 
parison. J] point out that generally (70) is much more 
exact than (71). The antenna current amplitude changes 
more rapidly in the beginning and there the series is 
least favorable. Observe also that in (70) and (71) the 
second negative power of 2,, resp ©, is lacking so that 


44E, Hallén, ‘Iterated sine and cosine integrals,” 


Trans. Roy 
Inst. Tech., Stockholm, no. 12, pp. 1-6; 1947, 
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the following terms are very small compared to the first 
one. Series corresponding to (70) for the outgoing wave 
and for reflected waves when the finite length of the 
antenna is taken into consideration have also earlier 
been derived (reference 8, page 1141 and page 1145; 
reference 9, eq. 55; reference 3, page 436). 

The series (70) is much more easily derived directly 
from the linearized integral equation than by series 
expansion of an exact solution (reference 3, pages 433~ 
436). The only attempt so far to evaluate for all z-values 
the exact expression (65b) for the wave amplitude of the 
reflected wave on a tube-shaped antenna has been series 
expansion, but as it involves the neglect of powers 
of af, this only leads to the already known approximate 
expression (70). 

It is, however, fortunate that at least the end value 
for 2 = 0 is exactly numerically known (the end ad- 
mittance) because when summing up the series of 
traveling waves on a finite antenna, the end admittances 
play an essential role and the total current is numerically 
sensitive just to this end value. So, even if at the moment 
the numerical evaluation of the exact reflected current 
wave along the whole antenna (from 65b) yet remains, 
the fact that we know W¢ (0) exactly is an essential result. 
Combined with (70) for those parts of the antenna which 
are not quite close to the end it gives a very exact 
knowledge of the antenna current everywhere. 

For the numerical evaluation of the exact expression 
for the antenna end admittance (50), (44) one has to 
transform the integral (44). This can be done in several 
ways. If we follow the same method as in (reference 3, 
page 431) for the linearized expression, we meet only two 
differences. They are 1) the poles of the transformed 
integral now form an infinite divergent series, and 2), 
the infinite integral, before compensation, diverges. 
Both difficulties can easily be overcome, and as this way 
of proceeding leads mainly to Bessel functions J,, 
Y,, of which very dense tables are available, we use this 
method. Otherwise, it is easy enough to avoid the series, 
but then the expression will retain functions J, and K, 
which are not so densely tabulated. 


I find (cf. reference 3, page 431) from (44): 


1 
2 SS SS 
$2 (8) —jr (a8) 


1 1+ Vi- wv? 
og 


1 
son {4 / Vesa ds 
"YS ouH (abu) H ( (apu) 
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oc u 
4 x log Vy — 1 du 
1 4H (abu). (apu) 
va8| du 


+ cl i ee | d — 
res 4 uv LuH \) (abu) Ho) (abu) 2 
B 


ng 2j FI 4 Sete log 2 |. (72) 
T 2 ap 


Except for the compensating last bracket term in the 
last infinite integral and the last terms outside the 
integrals, which all contain the factor af, one recognizes 
in this expression the corresponding linearized one 
(reference 3, eq. 35.58c) if as usual the Hankel functions 
are substituted by their first logarithmical part with 
suppressing of powers of a8. The integrals in (72) are 
easily converted into others with finite integrands and 
finite limits (cf. reference 3, eq. 35,.58d). The exact 
value of the end admittance of the tube shaped antenna 
(50) as it follows from (72) is shown in Fig. 1, together 
with the corresponding linearized (68), (69) and the series 
value as derived from (71). Fig. 2 gives a part of the 
exact curve in another scale. 

For the numerical evaluation of the end impedance 
of the inner current (63) we proceed in the same way, 
although the formula now differs because of the imagi- 
nary value of a in ; in (63). For summation of the series 
we need an asymptotic formula for 


n (Svs) 


for higher values of m. This asymptotic value is found 
to be 


J 2 2 ss Rop —jsgj st) 
(2 Vin? — aR) = Ves (cos j sin ; (73) 


m m™ 


J (a8) 
= Y ,(aB) 


2 if 1 u du \ 
mrJod+ Vi = uw?) H 0) (aBu)H > (apu)) 
&,, is defined by Jo(&m) = 0. 

The curves, shown in Fig. 1 and Fig. 2, speak for 
themselves. When the ratio of wave length to radius of 
cross section is more than 80-100 there is no visible 
difference whatsoever whether we use my old integral 
equation for thin antennas (the linearized equation) or 
the exact equation for tube-shaped antenna. 


CN? 


